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A  TEST  STRATEGY  FOR  HIGH  RESOLUTION  IMAGE  SCANNERS 


1.  Introduction 

The  scope  of  this  study  is  the  design  of  a  device  independent  test 
strategy  for  performance  tests  of  the  AIDS  (Advanced  ^mage  Digitizing 
^stem)  scanner.  Major  system  characteristics  to  be  tested  are  the  geo¬ 
metrical  accuracy,  linearity  of  gray  shade  response,  MTF,  and  resolution 
of  the  output.  The  study  is  of  a  purely  conceptual  type,  yet  most  of  the 
procedures  suggested  here  have  proven  their  practical  potential  elsewhere. 
The  major  performance  parameters  of  the  AIDS  are  described  in  AIDS  (1982). 
The  input  material  is  photographic  film.  The  AIDS  shall  be  capable  of 
scanning  over  a  maximum  format  of  9.5"  x  9.5".  The  minimum  scanned  image 
size  shall  be  1024  x  1024  pixels.  The  scanning  resolution  shall  be  100 
line  pairs  per  millimeter  over  a  format  of  6.6"  x  6.6"  at  a  MTF  level  of 
0.3.  As  a  secondary  requirement  the  AIDS  shall  scan  the  full  format  at 
larger  pixel  sizes,  corresponding  to  a  resolution- of  70,  50,  and  33  line 
pairs  per  millimeter,  at  a  MTF  level  of  0.3.  The  relative  geometric 
positional  accuracy  of  the  output  shall  be  one  pixel  mean  square  error 
over  the  entire  9.5"  x  9.5"  format. 

The  system  shall  provide  two  scanning  dynamic  ranges;  a  range  span¬ 
ning  1.25  density  units  and  a  range  spanning  2.50  density  units.  Pixels 
shall  be  generated  in  both  a  8-bit  and  a  11-bit  mode.  The  8-bit  mode  shall 
produce  data  having  255  equidensity  steps  with  one  step  repeatability. 
Linearity  shall  be  maintained  within  0.02  density  units  or  one  percent 
of  the  scanning  density  range,  whichever  is  greater.  The  11-bit  "inten¬ 
sity"  mode  shall  produce  1023  equi- intensity  steps  with  eight  steps  re¬ 
peatability.  The  linearity  shall  be  maintained  within  0.5  percent  of  the 
scanning  intensity  range. 

The  highest  average  scanning  rate  of  the  system  shall  be  not  less 
than  2  Mega  Pixels  per  second. 
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The  tests  in  this  study  are  designed  such  that  they  can  be  perfonned 
"in  situ",  i.e.  with  the  test  standards  placed  on  the  film  holder  and 
with  all  relevant  system  components  in  their  normal  position  of  operation. 

The  structure  of  this  report  distinguishes  clearly  between  tests  of 
the  geometric  and  the  radiometric  performance.  Both  parameters,  the  geo¬ 
metric  and  the  radiometric  accuracy,  however,  are  closely  related  in  prac¬ 
tice.  A  defect  in  the  radiometric  performance  of  an  electro-optical  ele¬ 
ment  might  well  influence  positional  measurements.  The  geometric  resolu¬ 
tion  is  dependent  on  the  radiometric  resolution  and  vice  versa.  The  MTF 
determination  with  multiple-bar  targets  across  scan  direction  might  be 
disturbed  by  tracking  errors  of  the  film  holder. 

Chapter  2  deals  with  the  geometric  performance.  Following  the  dis¬ 
cussion  of  some  possible  instrument  errors  a  number  of  test  standards  is 
suggested  for  use  in  the  test  procedures.  The  description  of  these  test 
procedures,  together  with  some  aspects  of  data  processing  and  data  analy¬ 
sis  concludes  Chapter  2. 

The  radiometric  performance  tests  are  depicted  in  Chapter  3.  Parti¬ 
cularly  touched  upon  are  MTF  measurement,  dynamic  range,  linearity  and 
noise,  response  uniformity,  flare  light,  temporal  stability,  and  coherency 
effects . 

Chapter  4  summarizes  the  test  procedures  in  table  format  and  indi¬ 
cates  options  and  priorities. 
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2.  Test  of  the  Geometrical  Performance 

The  Image  scanner  produces  a  digital  Image  of  the  scanned  analog 
film  Image.  The  geometric  performance  describes  the  degree  of  closeness 
(geometric  fidelity)  between  the  analog  and  the  digital  Image  (Input 
versus  output).  In  order  to  check  this  closeness,  test  patterns  of  known 
form  and  dimension  have  to  be  scanned.  A  comparison  between  the  Input 
pattern  and  the  output  pattern  shows  the  geometric  performance  of  the 
scanner.  The  major  problems  here  are  the  design  of  suitable  test  pat¬ 
terns  ("test  standards"),  the  development  of  test  procedures  and  data  pro¬ 
cessing  methods.  The  tests  of  the  geometrical  per  mance  can  be  formu¬ 
lated  such  that  they  are  resolution  independent.  ly  the  dimensions  of 
the  test  standards  must  be  selected  such  that  they  .low  for  sufficient 
responses  at  all  resolution  levels.  Critical  for  .-''cessful  geometrical 
test  Is  the  design  of  these  test  standards.  Test  '  ^ards  are  usually 
constructed  such  that  they  allow  for  the  detection  of  any  expected  Instru¬ 
ment  errors.  Since  the  constructional  details  of  the  AIDS  are  not  known 
at  this  point,  the  Instrument  errors  which  might  cause  geometrical  defor¬ 
mations  cannot  be  predicted  in  advance.  A  thorough  analysis  of  all  ex¬ 
pected  instrument  errors  is  not  feasible  at  this  stage.  The  test  stand¬ 
ards  are  selected  In  a  general  way  such  that  they  allow  for  the  detection 
of  possible  instrument  errors.  Some  possible  instrument  errors  are  de¬ 
scribed  in  the  following  chapter.  Most  of  these  errors  assume  a  linear 
array  to  be  used  for  scanning.  If  this  does  not  apply  the  correspondingly 
designed  tests  are  not  relevant. 


2.1  Some  possible  1nst»*ument  errors 

In  the  following  discussion  we  will  assume  that  the  scanner  consists 
of  a  linear  array  of  a  (variable)  number  e  (e  *  1,  .  .  n)  photosensi¬ 
tive  elements. 
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A  =  length  of  array 

a 


j linear  array  (e  elements) 


scan  direction  (center  line) 


Figure  1 :  Linear  array,  mounted  perpendicular  to  scan  direction 


2.1.1.  Dynamic  errors  in  straight  line  scanning 

Hoehle  (1983)  describes  three  major  types  of  errors  in  straight  line 
plotting  (Figure  2). 


wav  mess 
stepping 
tracking 


Figure  2:  Major  errors  in  straight  line  plotting 
(according  to  Hoehle  ;i983)) 

« 

Similar  errors  can  be  expected  in  a  general  scanning  process.  Waviness 
and  stepping  are  due  to  mechanical  oscillations  and  deformations  and  to 
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the  increment  size.  Stepping  should  only  appear  if  the  scanning  lines 
are  inclined  to  the  driver's  axes,  which  is  not  likely  to  be  the  case 
in  the  AIDS.  Waviness  and  tracking  errors  however  might  emerge  in  any 
scanning  system.  Tracking  errors  could  show  up  as  an  interval  error 
Ad  when  the  transition  from  one  line  to  another  is  made  (Figure  3  a)  or/ 
and  as  a  convergence/divergence  of  the  supposedly  parallel  scan  lines 
(Figure  3  b).  The  size  of  the  tracking  errors  Ad,  y  might  even  change 
from  one  line  pair  to  another. 
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a)  interval  error 


b)  convergence  error 


Figure  3:  Tracking  errors 


Besides  causing  geometrical  displacements,  waviness  and  tracking  errors 
might  produce  gaps  in  coverage  or  oversampling  in  and/or  across  scan  direc¬ 
tion. 


2.1.2.  Skewness  of  array  line 

Conceptually  the  linear  array  is  mounted  perpendicularly  to  the 
direction  of  scan  (Figure  3  a).  The  scan  interval  d  is  thus  equal  to 
the  length  of  the  array  A  skewness  of  the  array  relative  to  the 
scan  line  might  cause  gaps  in  the  scan  coverage  (Figure  4  b). 


zpc 


a)  correct  mount  b)  mounting  error  e,  caus¬ 
ing  gaps  across  scan 
direction 


c)  mounting  error  e, 
causing  oversampling 
in  scan  direction 


Figure  4;  Skewness  of  array  line 

wa  =  scanner  width  =  advancement  of  scanner 
zpc  ...  zones  of  partial  coverage 
zos  ...  zones  of  oversampling 


Certain  zones  (zpc)  are  only  covered  by  a  portion  of  an  array  element. 
The  width  of  the  zone  of  partial  coverage  is 
2AK  =  wg  sine  . 

The  interval  error  Ad  is 

Ad  =  tg  (1  -  cose)  . 

A  gap  in  coverage  results  if 


Ad  >  2AK  , 

or  ta  Cl  ■  cose)  >  wg  sine  , 

I  ^  ^  • 

For  a  small  e  we  get  the  condition 


e  > 


[ 
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A  loss  of  Information,  however,  might  occur  for  even  smaller  mounting 
errors  e,  because  along  zpc  the  size  of  the  effective  array  area  is 
reduced.  Oversampling  might  take  place  along  the  zones  zos  (compare 
Figure  4  c). 


2.1.3.  Length  of  array  line 

According  to  tests  conducted  at  Messerschmidt,  Boelkow,  Blohm  (0. 
Hofmann,  personal  communications)  linear  arrays  of  about  20  ran  length 
showed  approximately  ±10  pm  deviation  from  their  correct  length.  This 
deviation  can  be  corrected  during  the  scanner  calibration  process.  Re 
maining  deficiencies  in  the  length  of  the  array  line  would  cause  gaps 
or  overlap  in  the  output. 


2.1.4.  Linearity  of  array  line 

According  to  0.  Hofmann  (personal  coranuni cations) ,  the  misalignment 
of  CCD  array  elements  is  generally  negligible.  A  misalignment  would  pro¬ 
duce  "sheared"  images,  like  ti«  example  in  Figure  5. 


original  sheared  image 


Figure  5:  Sheared  image  caused  by  array  element  misalignment 


Misalignment  could  result  in  an  overall  deformation  of  the  line  (Figure  6 
or  could  show  up  as  a  sporadic  displacement  of  one  or  more  elements  (Fig¬ 
ure  6  b). 


a)  Overall  deformation 
Figure  6:  Types  of  misaligrenent 


b)  Displacement  of  a  set  of 
elements 


2.1,5.  Sampling  rate  errors  in  scan  direction 

If  the  sampling  rate  is  not  adjusted  to  the  physical  length  of  the 
pixel  in  scan  direction  and/or  if  any  constant  differences  between  these 
two  parameters  are  not  considered  in  the  definition  of  the  pixel  location 
a  geometrical  distortion  in  scan  direction  would  result.  In  the  simplest 
case  of  constant  difference  an  affine  deformation  would  ensue. 


2.1.6.  Other  system  errors 

Some  additional  errors  which  might  occur  in  the  system  are 

-  expansion  and  creep  of  the  optical  system  generated  by  the  heat  of 
the  light  source 

-  motion  of  film  caused  by  vibrations  at  high  scanning  speeds 

-  non-orthogonality  of  film  holder  to  the  optical  axis 

-  flatness  of  film  with  respect  to  holder 

-  relief  of  film/ emulsion. 
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These  errors  are  not  easily  separable,  but  like  the  others  they  are  part 
of  the  overall  error  budget. 


2.2.  Test  standards  and  test  patterns 

Test  standards  are  firmly  established  models  for  comparison  which 
serve  certain  test  purposes,  and  are  used  in  a  specific  way  in  well  proven 
test  procedures.  Test  patterns  imply  the  design  aspect  of  test  figures, 
and  are  used  as  such,  in  a  more  flexible  manner.  Since  there  are  no  es> 
tablished  test  procedures  available  for  our  test  task  we  will  not  dis¬ 
tinguish  between  the  terms  "standard"  and  "pattern".  The  selection  of  the 
test  figures  and  the  corresponding  test  procedures  should  be  done  such 
that  most  of  the  previously  mentioned  errors  (Chapter  2.1.)  can  be  checked. 
The  errors  listed  under  Chapter  2.1.6.  are  not  specifically  addressed. 


2.2.1.  Standard  test  grid 

This  is  the  most  popular  test  figure  for  geometrical  tests  of  x,  y- 
measurement  systems.  As  shown  in  Figure  7,  this  device  consists  of  a  num¬ 
ber  of  intersecting  parallel  lines.  The  grid  points  are  regularly  distri¬ 
buted  and  calibrated.  Wild,  Heerbrugg  offers  a  very  precise  glass  grid 
plate  with  a  coordinate  precision  of  the  calibrated  grid  points  of  0.5  urn. 
The  grid  points  are  distributed  over  an  area  of  24  x  24  cm2,  j^e  line 
width  is  constant.  The  grid  spacing  is  2  cm  in  both  directions.  It  allows 
for  a  check  of  the  global  positioning  accuracy  of  the  scanner.  For  details 
see  Appendix  A.l. 


Similar  to  the  standard  test  grid  this  device  consists  of  two  sets 
of  parallel  lines  intersecting  each  other.  However,  the  line  widths  are 
varying,  the  line  distances  are  smaller  and  the  overall  size  of  the  pat¬ 
tern  is  smaller. 


Figure  8:  Square  cross-line  graticule;  variable  line  width,  narrow  spacing 


This  device  is  readily  available  from  Heidenhain  in  nine  different  ver¬ 
sions.  Three  of  them  (nos.  14,  15,  16  for  the  different  scanner  resolu 
tion  levels)  could  be  considered  here.  See  Appendix  A. 2. 


2.2.3.  Dot  and  line  comparison 

Appendix  A. 3.  shows  the  dot  and  line  comparison  chart  no.  20  from 
Heidenhain.  The  widths  of  the  lines  and  the  diameters  of  the  dots  range 
from  10  urn  to  1000  um.  The  spacing  of  the  lines/dots  of  different  widths 
is  2  mm.  This  chart  can  be  used  to  check  the  resolution  of  individual 
fundamental  patterns  like  single  lines  and  dots.  It  can  also  be  applied 
to  check  the  local  geometric  performance  of  the  scanner  with  respect  to 
the  rendering  of  these  specific  objects. 


2.2.4.  Parallel  line  test  chart 

The  line  test  chart  no.  32  from  Heidenhain  is  shown  in  Appendix  A. 4. 
This  chart  with  its  variable  line  pitch  can  be  used  to  check  the  length 
of  the  sensor  array. 


2.2.5.  Radial  circle  test  chart 

This  test  chart  is  shown  in  Figure  9.  It  is  not  readily  available, 
but  it  could  possibly  be  manufactured  by  Heidenhain.  This  chart  has  the 
very  desirable  property  of  isotropy.  The  location  of  the  center  of  the 
circles  in  the  scanner  output  gives  an  indication  of  the  local  position¬ 
ing  accuracy.  In  addition  the  rendering  of  another  fundamental  object 
pattern,  the  circle,  can  be  checked. 
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Dimensions  of  radial  circle  test  chart  figures: 


Distances  of  circles: 

Si  »  0.18  nw 
$2  *  0.36  tim 
S3  *  0.72  mm 
Si*  *  1.44  mm 


Diameters  of  circles: 

Di  =  35  ym 
D2  =  70  ym 
D3  *  140  ym 
D4  =  280  ym 
Ds  =  560  ym 


Line  widths: 

di  =  5  ym 
d2  *  10  ym 
ds  *  20  ym 
d4  *  40  ym 
ds  =  80  ym 


Lengths  of  cross  bars 

tj  =  20  ym 
I2  -  40  ym 
li  =  80  ym 
2,4  *  160  ym 
Is  -  320  ym 
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2.3.  Test  procedures,  data  processing,  and  data  analysis 

According  to  the  purchase  description  of  the  AIDS  (AIDS,  1982)  a 
display  monitor  and  a  cursor  will  be  available  to  allow  the  operator  to 
perform  measurements.  The  measurement  resolution  shall  be  one-third 
pixel,  which  is  well  below  the  required  one  pixel  positioning  accuracy. 
Since  the  scanner  output  data  is  in  digital  form  the  measurements  could 
either  be  performed  in  the  digital  domain  or  softcopies  could  be  measured 
using  the  monitor/cursor  system.  The  latter  option  is  preferred  for  the 
geometrical  tests  because  it  does  not  require  the  extensive  software  which 
would  be  necessary  for  the  interpretation  and  measurement  of  digital  data. 
Any  operator  related  subjectivity  in  the  measurement  process  is  part  of 
the  system  and  would  occur  in  the  measurement  of  real  photographs  as  well. 
For  the  test  purpose  at  hand  an  operator-driven  measurement  system  is  opti¬ 
mal,  it  is  flexible,  fast,  and  inexpensive.  For  the  radiometric  tests  it 
is  crucial  to  operate  directly  in  the  digital  domain.  Any  intermediate 
analog  display  of  the  digitized  values  would  cause  an  additional  modula¬ 
tion  and  would  require  cascading. 

The  test  figures  suggested  in  Chapter  2.2.  are  gratings  and  grati¬ 
cules  on  glass  substrate.  Since  the  AIDS  can  only,  scan  film  based  ob¬ 
jects  the  test  figures  have  to  be  copied  on  film.  It  is  suggested  to  use 
Kodak  High  Resolution  Duplicating  Film  (Ester  Thick  Base)  SO-187  (see 
Appendix  B.I.).  The  film  based  test  figures  have  to  be  measured  on  a 
precision  stereo-comparator  or  analytical  plotter,  stored  and  used  under 
constant  humidity  and  temperature  conditions.  A  possible  distortion  of 
the  test  patterns,  caused  by  a  deformation  of  the  film  base/emulsion, 
should  be  monitored  through  occasional  remeasurements. 
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2.3.1.  Test  of  the  global  geometric  perfomance 

This  test  uses  the  standard  test  grid  described  under  2.2.1.  in 
order  to  determine  the  overall  (global)  geometric  performance.  The  scan 
ner  output  of  the  13  x  13  =  169  crosses  of  the  test  grid  is  to  be  observ 
ed  in  two  independent  sets  of  measurements.  It  is  suggested  that  the 
measurements  be  modelled  by  a  bivariate,  orthogonal  polynomial.  Bivar¬ 
iate  orthogonal  polynomials  do  have  the  following  advantages  compared  to 
conventional  polynomials: 

-  functional ;  bivariate  polynomials  compensate  for  the  total 
deformation  in  all  points  of  a  regular  grid  if  the  order 
of  the  polynomial  is  adjusted  to  the  number  of  grid  points 

-  numerical ;  orthogonal  polynomials  produce  optimal  system 
stability;  variances  and  covariances  are  minimal  (Rao  1973, 
p.  235  ff) 

-  statistical;  orthogonal  polynomials  lead  to  a  higher  power 
(type  II  error)  in  significance  testing  (Graybill  1961,  p. 

140) 
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Figure  10:  5  x  5  =  25  grid  point  distribution 

A  bivariate  polynomial,  which  is  orthogonal  with  respect  to  a  5  x  5  =  25 
regular  grid  point  distribution  (Figure  10)  can  be  formulated  as 

Ax  »  y^ As7, 

Ay  ■  y^BjX 


(1) 
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with 


*  (1  y^  (y^-2b^)  y(y^-  y*{y^-  +  Hb"* 

“^  =  (1  x"  (x2-2b2)  x(x"- -^^2)  ^2(x2.  Ikzx  .  72.  u 


X'  = 


and 


an  a^  2  •  ♦  •  ai  5 

bjj  bj2  ...  bjj 

•  •  • 

•  •  • 

• 

A,  “ 

•  •  • 

•  •  • 

.  B3  = 

•  •  • 

•  •  • 

.asi  a52  ...  a55_ 

.bsi  ^5  2  ...  bss^ 

With  the  substitutions 

k  »  x^-  2b*,  i  =  y*-  2b*,  p  »  x*-  Q  *  y*- 

r  -  xHx‘-  +  3|b\  s  • 


we  obtain 
Ax  * 


*  an+  aijX  +  aj^y  +  ajjXy  +  a^jk  +  aj^il  +  a^jyk  +  ajjxA  + 
+  ajjkA  +  aiuxp  +  a^yq  +  aj^xyp  +  a^jxyq  +  aj^xtp  + 

+  a^sykq  +  aH^xypq  +  aisr  +  +  a^jyr  +  ajjxs  + 

+  a35  r  +  ajjks  +  a^jyqr  +  aj^xps  +  ajjTS  ; 

^y  =  ‘>11+  ‘>i2>‘  +  b^y  +  bjjxy  +  b„k 

+  bjjkt  +  bj^xp  +  b^yq  +  bj^xyp  +  b^jxyq  +  bj^x^p  + 

+  b43ykq  +  b^^xypq  +  buT  +  bj^s  +  b^jyr  +  bj^xs  + 

+  b35  r  +  bjjks  +  b^jyqr  +  b^^xps  +  b^^rs 


(2) 


(3) 


(4) 


[ 


(5) 


If  the  number  of  grid  points  is  increased  beyond  25  the  property  of  ortho¬ 
gonality  is  lost,  but  even  in  this  case  the  suggested  polynomial  is  closer 
to  orthogonality  than  a  conventional  polynomial. 


If  n  grid  points  are  used  for  polynomial  coefficient  determination,  the 
average  values  of  the  observed  grid  point  coordinates  are  computed  as 


I 
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Systan  C8)  consists  actually  of  two  independent  estimation  models,  which 
can  be  adjusted  separately,  as  long  as  the  joint  weight  matrix  P  has  the 
blockdiagonal  structure  as  indicated.  In  fact  we  are  going  a  step  fur¬ 
ther  in  assuming  that  P^*  Py  have  diagonal  structures  with  the  diagonal 
elements  derived  from  equation  (7).  The  least  squares  approach  leads 
thus  to  two  sets  of  normal  equations 


with 


N  p 
x'^x 


NyPy 


=  (a'p  a  ) 

XXX 

=  (AjP/,) 

^x^x^^ 

=  AjpyAy 


(9) 


Analysis  of  results: 

The  adjustment  procedure  will  deliver  the  following  values. 
Estimated  polynomial  parameters: 

Px  =  ^x'  ‘-x  ’ 

Py  -  Ly  . 


(10) 


Residuals : 


Vy  •  AyPy  -  iy  . 


(U) 


i 
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Variances  of  unit  weight: 

.V2  ,  VxPxVx 
%x  ~  n  -  njj  * 


wT  p  w 

y  y  y 

n  -  n.. 


nj(,  Oy  . . .  number  of  polynomial  para¬ 
meters  in  Pj(,  Py 


(12) 


Covariance  matrices  of  the  estimated  polynomial  parameters: 


N'' 
ox  X 


(13) 


Significance  testing  of  polynomial  parameters: 
The  test  criteria 


T 


X 


n 


T 

y 


/<s  T  *  1  ^ 
P  Z  ^  p 

n.. 


(14) 


can  be  used  to  test  whether  the  whole  set  of  polynomial  parameters  is 
significant  or  not.  Under  the  null -hypotheses 


(15) 
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the  test  criteria  T 


X’ 


Ty  follow  the  central  Fisher  distributions 


'  F  "-"x) 

Ty  -  F  C"y.  n-ny) 

We  reject  if 


(16) 


T  >  C  =  F  (l-a,  n  ,  n-n  ),  a  ...  type  I  error  size  of  the 

global  tests 

Ty  >  Cy  =  F  (l-a.  Hy,  n-ny).  (17) 


If  any  one  of  the  global  null -hypotheses  (15)  is  rejected,  the  signifi¬ 
cance  of  the  associated  individual  polynomial  parameters  should  be  tested. 
The  probability  levels  of  these  individual  tests  are  related  to  those  of 
the  general  test.  If  the  individual  polynomial  parameters  would  be  inde¬ 
pendent  of  each  other,  the  type  I  error  size  a  of  the  global  tests  would 
be  related  to  the  type  I  error  sizes  a^.  Uy  of  the  individual  tests  ac¬ 
cording  to: 


(l-a)  =  (1-a^)'’’^ 
(l-a)  =  (l-ay)"y 


(18) 


Strict  orthogonality  among  the  polynomial  parameters  is  not  valid  any 
more  if  more  than  the  necessary  number  of  grid  points  (25  points  in  this 
example)  is  used.  A  solution  to  this  problem  is  to  orthogonal ize  the 
vector  of  polynomial  parameters  a-posteriori  and  to  test  the  orthogonal - 
ized  components.  This  can  be  achieved  practically  without  any  additional 
computational  amount,  if  the  solution  of  (9)  is  obtained  by  triangular 
factorization  of  Nx,  Ny.  The  intention  is  to  transform  the  original 
solution  vectors  Pj^,  Py  (distributions:  Px^N  (Pj^,  ip  ), 

Py-N  (Py,  Spy  ))  with  the  transformations 
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(19) 


into  orthogonal  components  fx,  fy  with  the  distributions  ’ 

fy~N(FyPy,  aQyQfy).  It  can  be  shown  that  F^  can  be  chosen  such  that 
Qfjj  =  Qfy  *  I,  so  that  the  components  of  fx*  fy  are  independent  on  each 
other.  If  the  normal  equations  matrices  Nx*  Ny  are  Choi esky-factori zed 
with 

Ny  ®  =  C  cj  ,  (20) 

X  XX  y  y  y 

the  system  (9)  takes  the  form 


or 


C  C^p  »  L  . 

X  x^x  X  * 


C  C^p  *  L 
y  y^y  ‘'y 


T>'  -1 

C  p  =  C  'l 

X^X  X  X  ’ 

cX-c;-L,  • 


(21) 


(.22) 


It  can  be  shown  (Finn  1974,  p.  138)  that  the  orthogonal  components  of 


A  A 


Px'  Py 


^“CLx  • 


f  =  C  L 
y  y  y  ’ 


(23) 


which  are  the  factorized  "right  hand  sides"  of  the  normal  equations  (9). 
These  vectors  are  generated  during  the  factorization  procedure.  All  that 
is  needed  is  to  test  their  components. 
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Under  the  null -hypotheses 


: 

f 

=  0 

ox 

iJ  : 

f 

=  0 

oy 

the  corresponding  test  criteria 

.  o~  ’ 

1  Ox 


=  . 

oy 


(24) 


(25) 


are  Student-distributed  with 

tx  t(n-nj^)  , 


If 


ty^  -t(n-ny)  . 


*  '^x  ’  "-"x> 


*y.  *  ‘^y  '  "'"y' 

j 


(26) 


(27) 


the  null -hypotheses  (24)  are  accepted. 

In  this  case  the  corresponding  components  ?xi*  ^yj  set  to  zero 
{^x^  “  0,  ^yj  =  0).  Finally  the  original  polynomial  parameters  are 
puted  by  backsubstitution  from 


com- 


C  p  =  f  , 

x^x  X 


C  p  *  f  , 
y^y  y 


(28) 
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whereby  f',  f'  are  the  "cleaned"  orthogonal  parameter  vectors.  All  other 
X  y 

aspects  of  the  data  processing  procedure  are,  for  the  most  part,  standard. 
Blunder  detection  can  be  performed  by  computing  the  residual  vectors  v^, 

Vy  and  applying  Baarda's  data-snooping  technique  (compare  e.g.  Gruen,  1978). 
It  may  be  of  interest  to  investigate  the  type  of  deformation  involved  if 
only  a  three-parameter  transformation  is  applied,  i.e,  two  translation 
components  and  one  rotation. 

To  achieve  this,  use  the  transformation: 


AX  =  y'^Aj^x 

Ay  =  y^Bg^x 


(29) 


with 


(30) 


The  adjustment  is  performed  as  previously  outlined,  except  that  the  test¬ 
ing  of  parameters  is  not  necessary  in  this  case.  Here  the  pattern  gener¬ 
ated  by  the  residuals  is  of  particular  interest.  Tt  is  suggested  that 
the  residuals  v^,  Vy  be  plotted.  A  visual  inspection  of  this  plot  will 
reveal  any  deformations.  A  scale  factor  in  the  data  can  be  compensated  by 
a  similarity  transformation.  Two  different  scale  factors  along  the  co¬ 
ordinate  axes  and  a  shearing  effect  can  be  accommodated  by  an  affine  trans¬ 
formation.  The  formulae  for  these  transformations  are 
similarity  transformation: 


Ax 


y^AgX 


Ay 


(31) 
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with 


■  an 

ai2l 

^21 

L  ®21 

- 1 

O 

Bs  - 

_ 1 

0  . 

affine  transformation: 
Ax  *  y^A  X 

d 


with 


Ay  =  y^BgX 


ail 

aial 

1 _ 

^1  2 

321 

- 1 

O 

• 

Ba  “ 

r>4 

_ _ 1 

0 

C32) 


(33) 


(34) 


2.3.2.  Skewness  and  linearity  of  the  array  line 

Both  errors  can  be  checked  by  scanning  one  or  nore  lines  which  are 
inclined  to  the  scan  direction.  The  standard  test  grid  and/or  the  square 
cross-line  graticule  can  be  used  for  this  purpose.  The  test  figure  should 
be  mounted  at  the  scanner  stage  such  that  one  group  of  lines  is  approxi¬ 
mately  perpendicular  to  the  scan  direction.  Visual  inspection  of  the 
output  for  these  lines  would  indicate  a  defect.  The  skewness  of  the  ar¬ 
ray  line  would  show  up  in  a  sawtooth  shaped  output  line  (Figure  11). 


scan  direction 


correct  output 


distorted  output 


Figure  11:  Effect  of  skewness  of  array  line  on  output  of  test  line  per¬ 
pendicular  to  scan  direction 
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A  non-linearity  of  the  array  line  would  result  in  a  distortion  of  the 
test  lines.  The  pattern  of  distortion  would  depend  on  the  type  of  non 
linearity.  It  would  show  up  as  a  regular  pattern  with  a  spacing  equal 
to  the  array  length  Csimilar  to  the  skewness  pattern). 


2.3.3.  Interval  tracking  error  and  length  error  of  the  array  line 

A  constant  interval  tracking  error  CChapter  2.1.1.)  and  a  length 
error  (Chapter  2.1.3.)  of  the  array  line  can  be  separated  from  each 
other  if  the  output  of  one  single  array  line  is  extracted  and  analyzed. 
In  order  to  check  the  length  of  the  array  line,  the  use  of  Heidenhain's 
test  chart  no.  32  (Appendix  A. 4.)  is  suggested.  This  chart  should  be 
scanned  with  the  lines  parallel  to  the  scan  direction.  The  difference 
between  input  and  output  line  spacing  gives  directly  an  indication  of 
an  array  line  length  error.  The  film  based  line  test  chart  should  be 
measured  on  a  precision  comparator  or  analytical  plotter  prior  to  the 
scanning  process  in  order  to  verify  the  line  widths  and  distances. 

For  the  check  of  the  constant  interval  tracking  error  the  standard 
test  grid  and/or  the  square  cross-line  graticule  can  be  used.  The  test 
figure  should  be  set  up  such  that  one  group  of  lines  is  parallel  to  the 
scan  direction.  The  comparison  between  input  and  output  distances  be¬ 
tween  these  parallel  lines  show  the  accumulated  effect  of  this  tracking 
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Figure  12:  Determination  of  the  constant  interval  tracking  error 

The  non-constant  portion  of  this  tracking  error  can  Jje  checked  if  sets 
of  two  adjacent  scans  are  analyzed.  This  can  be  a  very  tedious  proce¬ 
dure  if  many  scan  lines  are  tested.  The  line  test  chart  no.  32  (Appendix 
A. 4.)  could  be  used  for  this  purpose.  If  not  removed,  the  error  would 
show  up  as  part  of  the  total  error  budget  as  computed  under  test  2.3.1. 


2.3.4.  Convergence/divergence  tracking  error 

By  monitoring  the  parallelism  of  the  output  of  two  parallel  input 
lines  this  type  of  error  can  be  analyzed.  The  input  lines  should  be 
parallel  to  the  scan  direction.  The  test  grid  lines  (Appendix  A.l.)  could 
serve  this  purpose.  Again  the  overall  tracking  error  might  be  split  up 
into  a  constant  portion  and  an  irregular  part.  The  constant  portion  can 
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be  determined  by  computing  the  average  y-value  (7)  for  each  set  of  two 
adjacent  lines  (see  Figure  13). 

scan  direction 


adjacent  grid  lines 


Figure  13:  Determination  of  the  constant  convergence/divergence  track¬ 
ing  error 

For  the  monitoring  of  the  non-constant  portion  we  face  a  similar  problan 
as  under  2.3.3.  Each  pair  of  two  adjacent  lines  must  be  investigated. 

A  possible  error  of  this  type  would  show  up  in  test  2.3.1. 


2.3.5.  Sampling  rate  error  in  scan  direction 

This  error  causes  a  scale  change  in  the  direction  of  scan  and  can 
be  detected  by  checking  the  output  distances  between  a  set  of  parallel 
lines  running  perpendicular  to  the  scan  direction.  The  standard  test 
grid  lines  can  be  used  for  this  purpose.  In  addition,  the  global  test 
2.3.1.  delivers  the  scale  factor  baj  in  y-(scan)  direction.  A  signifi¬ 
cant  deviation  of  bjjfrom  the  value  1  would  indicate  a  possible  sampling 
rate  error. 
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2.3.6.  Stepping  and  waviness 

These  errors  can  be  detected  simply  by  checking  the  straightness 
of  a  standard  test  grid  line  in  scan  direction.  A  qualitative  check 
can  be  performed  visually.  Measurements  should  be  made  only  if  these 
errors  have  to  be  quantized. 


2.3.7.  Local  performance  and  rendering  of  fundamental  patterns 

In  addition  to  a  high  geometrical  positioning  accuracy  it  is  of 
vital  importance  for  a  precision  scanner  to  render  certain  features  ac- 
cu 'ately  and  completely.  The  most  essential  geometrical  patterns,  which 
form  a  basis  for  other  more  complicated  features  imaged  on  aerial  photo¬ 
graphs,  are  dot,  circle,  and  bar.  It  is  suggested  that  the  output  of  the 
AIDS  be  checked  with  respect  to  the  rendering  of  these  patterns. 

The  dot  line  comparison  chart  (Appendix  A. 3.)  and  the  radial  circle  test 
chart  (Figure  9)  can  be  used  as  test  standards..  The  radial  circle 
test  chart  can  also  serve  as  a  standard  to  check  the  local  geometrical 
positioning  accuracy.  The  standard  test  grid  with  its  grid  width  of  2  cm 
can  only  control  the  global  accuracy.  The  radial  circle  test  chart 
with  its  33  circles  of  different  diameters  and  line  widths,  and  with  the 
center  crosses  spread  point-symmetrically  over  a  6  x  6  rnn^  area,  can 
control  the  local  positioning  accuracy.  If  the  system  works  correctly, 
a  3-parameter  transformation  (equations  (29),  (30))  should  be  sufficient 
to  connect  the  output  coordinates  of  the  crosses  to  the  given,  local 
values.  Any  additional  significant  polynomial  terms  indicate  a  malfunc¬ 
tioning. 
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3.  Test  of  the  radiometric  performance 
3.1.  Introduction 

The  three  basic  image  quality  parameters  that  are  used  to  describe 
the  performance  of  a  high  speed  film  scanner-digitizer  are  MTF,  linearity 
of  response  and  noise  level  over  the  operating  or  dynamic  range  of  che 
system.  An  additional  parameter  is  the  amount  of  flare  or  stray  light 
present  in  the  system.  This  was  not  specified  in  the  AIDS  purchase  de¬ 
scription  but  will  be  discussed  later  in  this  report.  Three  quantities 
that  are  seldom,  if  ever,  specified  are  response  uniformity,  temporal 
stability  and  the  influence  of  coherency  effects  on  performance.  These 
will  also  be  discussed  later  in  this  report. 

When  considering  MTF  and  response  linearity  one  must  be  careful  to 
evaluate  the  system  at  points  throughout  its  scanning  speed  range  as  a 
function  of  target  spatial  frequency  and  for  all  outputs.  (There  will 
be  only  one  output  for  a  conventional  microdensitometer  type  scanner, 
but  there  could  be  2048  or  more  for  a  CCD  system) .  Of  particular  concern 
here  is  the  frequency  response  of  the  detector  electronics,  in  particular 
the  analog  to  digital  (A/0)  converter.  Problems  have  been  encountered 
with  a  scene  dependent  lag  in  response  in  the  along-scan  direction  due 
to  the  Butterworth  filter  in  the  output  of  NASA's  multi -spectral  scanner 
system;  four  detectors  in  band  5  of  the  thematic  mapper  (TM)  show  an  over¬ 
shoot  response  when  crossing  a  white-dark  boundary,  (Kieffer  et  al.,  1983), 
and  the  PDS  microdensitometer  (Slater,  1980,  pp.  388-394)  exhibits  a 
significant  reduction  in  modulation  output  at  scan  speeds  above  1  mm/s 
when  scanning  a  316  cy/mm  target  with  a  1  ym  scanning  aperture,  see  Fig¬ 
ure  19.  The  bin  sizes  (output  level  steps  of  the  analog  to  digital  con¬ 
verter)  in  the  TM  have  been  found  to  be  unequal,  which  results  in  the 
histograms  of  uncorrected  data  from  a  random  scene  having  a  comb-like 
appearance  rather  than  being  smooth  in  shape.  This  problem  reduces  the 
quantization  resolution  from  the  advertized  8-bit  to  effectively  7-bit 
CKieffer  et  al.,  1983). 
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Scans  to  determine  the  effect  of  scan  speed  on  recorded  modu¬ 
lation.  Conditions:  lOX  objectives;  white  light;  200X  magni 
fication;  "A"  apertures;  1-um  scanning  aperture.  Target 
frequency  316  cycles/mm 
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Although  the  reasons  for  some  of  these  problems  have  not  yet  been  fully 
determined,  they  relate,  in  the  TM  case,  to  the  difficulty  of  sampling 
and  digitizing  some  4  x  10®  point  values  per  second.  This  is  of  the  same 
order  as  the  AIDS  and  two  orders  of  magnitude  higher  than  the  PDS  micro¬ 
densitometer.  A  second  concern  is  that  the  response  of  the  system  be 
determined  at  various  locations  over  the  entire  film  format  in  order  to 
check  focus.  The  focus  problan  divides  into  four  parts:  First,  the 
orthogonality  of  the  film  holder  to  the  optical  axis,  second,  the  flat¬ 
ness  of  the  film  with  respect  to  the  holder,  third,  the  parallelism  of 
motion  of  the  film  holder,  and  fourth  the  quality  of  focus  along  the  length 
of  the  scanning  aperture.  With  an  autofocus  device,  the  first  three  may 
be  insignificant  but  the  sensitivity  and  frequency  response  of  the  auto- 
focusing  need  to  be  determined. 


3.2.  MTF  measurement 

There  are  two  common  ways  to  determine  system  MTF;  the  first  is  by 
reference  to  a  square  bar  target  array,  the  second  is  by  edge  gradient 
analysis.  The  former  is  simpler  in  terms  of  post-scan  data  processing 
but  is  more  expensive  because  of  the  high  quality  resolution  targets  re¬ 
quired,  Both  methods  will  be  described  briefly  in  this  section. 

The  bar  target  MTF  can  most  conveniently  be  determined  by  use  of  a 
mul^ibar  target.  Typically  these  are  high  or  low  contrast  targets  on 
film  or  glass  that  consist  of  sets  of  15  light  bars  on  a  dark  background 
as  shown  in  Figure  15. 
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Figure  15:  Example  of  15-bar  target 

The  choice  of  a  15-bar  target  is  made  for  two  reasons;  first,  it 
provides  about  10  cycles,  which  are  not  distorted  because  of  being  near 
the  ends  of  the  pattern,  to  be  averaged  to  get  an  accurate  determination 
of  modulation.  Second,  when  scanned  slowly,  the  focus  can  be  easily  ad¬ 
justed  to  provide  maximum  modulation.  This  electrical  focusing  is  more 
sensitive,  repeatable  and  provides  a  higher  modulation  than  visual  focus¬ 
ing.  Three-bar  targets  or  continuously  varying  frequency  targets  do  not 
have  these  advantages.  Fifteen-bar  targets,  similar  to  those  shown  in 
Figure  15,  can  be  produced  on  film.  A  set  of  standard  3-bar  targets  on 
9.5  inch  wide  film  in  an  arrangement  suitable  for  AWAR  (Area  Weighted 
Average  Resolution)  testing  has  recently  been  sold  to  DMA  by  Itek  Corp. 
according  to  Ron  Carbonier  of  Itek  (phone  (617)  276-2551).  He  estimated 
that  a  similar  arrangement  of  15-bar  targets  covering  the  spatial  fre¬ 
quency  range  1-300  cy/mm  would  cost  of  the  order  of  $2000. 

In  the  particular  example  shown  in  Figure  15,  in  all  but  one  case  there 
are  frequencies  close  to  the  specified  values  of  33,  50,  70  and  100  cy/mm. 
The  .MTF  is  found  by  scanning  the  target  at  the  various  frequencies  up  to 
as  much  as  twice  the  frequency  of  interest.  The  modulation  at  each  fre¬ 
quency  is  determined  by  the  standard  equation  and  the  MTF  curve  is  plotted 
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technique  has  been  used  extensively  in  the  reconnaissance  conriunity  (e.g. 
in  2131  to  determine  overall  in-flight  system  MTF.  Thorough  descriptions 
can  be  found  in  Thomas  {.1973),  pp.  966-972  and  Dainty,  Shaw  (1974),  pp. 
244-246. 


Figure  17:  Derivation  of  the  MTF  from  the  edge  response  curve  (after 
Dainty,  Shaw,  1974,  p.  245) 


The  first  step  in  the  procedure  is  to  produce  an  extremely  sharp 
edge  between  two  large  uniform  areas  of  different  densities.  Usually 
the  edge  of  a  razor  blade  is  contact  printed  onto  an  extremely  high  re¬ 
solving  power  film  such  as  Eastman  Kodak  649F,  brush  development  is  used 
to  minimize  adjacency  effects.  To  calibrate  the  test  edge  function,  a 
scan  can  be  made  with  a  1  ym  slit  and  a  1  urn  sampling  interval  using  a 
high  performance  microdensitometer.  A  typical  result  is  shown  in  Figure 
18  of  the  spatial  frequency  of  the  edge  as  determined  by  an  TGA  computer 
program  (written  by  R.A.  Schowengerdt  at  the  University  of  Arizona). 
Corrections  of  this  spectrum  should  be  made  for  the  width  of  the  scan¬ 
ning  slit  using  the  sine  function  for  the  slit. 

In  the  second  step,  the  scanner/digitizer  is  used  to  scan  the  edge 
several  times  in  both  directions  at  different  speeds  to  determine  possible 
hysteresis  effects  (where  the  response  is  modified  by  the  proceeding 
signal  level)  or  frequency  response  deficiencies  in  the  detector(s)  and 
electronics.  For  this  measurement,  the  scan  direction  should  be  accurately 
perpendicular  to  the  straight  edge. 
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Figure  18:  Calibration  edge  spectrum  multiplied  by  microdensitometer 
OTF  for  0.001  by  0.100-mm  aperture 


The  third  step  is  manually  to  draw  a  smooth  curve  through  the  edge 
trace  or  use  a  computer  program  to  best-fit  polynomials  over  the  toe  and 
shoulder  of  the  curve.  If  the  latter  program  written  by  Schowengerdt  is 
used  it  will  continue  to  plot  the  required  square  wave  MTF,  at  whatever 
spatial  frequency  intervals  are  requested,  by  differentiating  the  smoothed 
edge  trace  and  then  taking  the  Fourier  transform  of  the  result.  A  manual 
approach  is  described  in  Thomas  (1973),  but  for  the  purpose  of  evaluating 
AIDS,  a  computerized  approach  is  recommended.  Schowengerdt 's  program  or 
software  at  213  should  be  available. 

Finally,  the  frequency  spectrum  of  the  edge  should  be  accounted  for 
by  dividing  the  MTF  obtained  in  step  3  by  that  spectrum. 

In  the  data  analysis  step,  if  a  CCD  linear  array  is  used,  care  should 
be  taken  to  examine  the  output  from  all  detectors.  There  can  be  signi¬ 
ficant  differences  from  detector  to  detector  as  well  as  slowly  varying 
changes  in  sensitivity  across  the  array.  Both  of  these  can  generally  be 


corrected  to  about  the  0.1%  level.  In  addition,  the  effect  of  scan  speed 
and  sampling  frequency  on  MTF  for  each  detector  should  be  determined. 

Even  if  the  EGA  method  is  not  used  for  MTF  determination,  the  edge  re¬ 
sponse  should  be  used  for  carefully  checking  response  uniformity  as  a 
function  of  scan  speed,  sampling  frequency  and  contrast. 

MTF  curves  obtained  by  the  EGA  method  for  a  PDS  microdensitometer 
are  shown  for  various  scan  speeds,  two  scan  directions  and  for  two  adjust 
ments  in  Figures  19,  20  and  21. 


Spatial  frequency  (cycles/mm) 


Figure  19:  Microdensitometer  OTF  for  different  scan  speeds  for  a  0.040-niii 
aperture 
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figure  20;  Microdensitometer  OTF  for  opposite  scan  directions  over  edg 
for  0.010-mm  aperture  (scan  speed  =0.3  im/s) 


Figure  21:  Microdensitometer  OTF  reproducibility  for  two  independent 
alignments  by  same  operator  for  0.010-mm  aperture 
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3.3.  Dynamic  range,  linearity  and  noise 

The  purpose  of  the  following  test  is  to  check  the  radiometric  re¬ 
sponse  of  the  AIDS  while  varying  parameters  such  as  dectector  gain,  a- 
perture  size  and  calibration  neutral  density  filters  (if  included). 

Two  types  of  targets  are  available  for  this  test.  One  is  a  grey 
step  tablet  available  from  Eastman  Kodak  that  varies  from  clear  to  a  den¬ 
sity  of  3.0  in  about  20  steps.  The  other  is  a  continuous  grey  scale  that 
can  be  fabricated  by  grinding  and  polishing  a  wedge  from  a  thick  glass 
neutral  density  filter  and  then  contact  printing  the  wedge  onto  film. 

Some  results  from  tests  on  a  PDS  microdensitometer  are  shown  in 
Figure  22.  Note  that  this  instrument  uses  a  photomultiplier  tube  that 
has  a  larger  dynamic  range  than  CCD  arrays  and  the  non-linearity  effects 
are  only  observable  at  very  high  and  low  light  levels,  when  the  applied 
PMT  voltages  were  low  and  high  respectively.  For  example,  for  the  very 
large  40  x  400  ym  slit  with  no  ND  filter  (the  fourth  line  down)  the  PMT 
voltage  had  to  be  set  very  low  so  as  not  to  exceed  the  maximum  signal 
out.  This  gave  rise  to  a  non-linear  response  over  the  0-1.0  target  den¬ 
sity  range.  At  the  other  extreme,  the  lowest  fan  of  three  lines  covers 
the  density  range  5  to  7.  Here  the  PMT  voltage  had  to  be  set  at  its 
highest  permissible  level  and  not  only  was  a  non-linear  response  observed 
but  also  the  effect  of  noise  in  the  width  of  the  fan.  AIDS  will  not  have 
to  operate  over  a  range  of  more  than  0  to  2.5  density  units  so  we  might 
anticipate  that  non-linearity,  noise  and  dynamic  range  would  not  be  a 
problem.  We  should  note,  however,  that  the  dynamic  range  of  CCDs  may  be 
less  than  1000:1  and  that  at  high  speeds  the  signal  to  noise  ratio  will 
be  reduced  because  of  the  decreased  integration  time.  Again  these  tests 
should  be  made  as  a  function  of  scan  speed  and  sampling  rate. 

The  dynamic  range  and  linearity  test  can  be  conducted  with  either 
the  step  tablet  or  the  continuous  wedge.  Either  should  be  checked  with 
a  PDS  microdensitometer  using  a  large  aperture  (say  40  x  200  ym) .  The 
linearity  of  the  PDS  can  easily  be  checked  by  scanning  a  straight  edge 
across  the  large  aperture. 


40 


The  use  of  the  readily  available  Eastman  Kodak  20-step  tablet  for 
the  linearity  test  at  the  0.5%  level  may  be  marginal.  A  tablet  contain¬ 
ing  a  larger  number  of  steps  can  be  readily  fabricated  in  a  photographic 
laboratory. 

To  determine  the  noise  level  of  the  system,  and  to  avoid  spatial 
variations  in  the  filter  being  misinterpreted  as  noise,  a  2.5  neutral 
density  filter  should  be  placed  in  the  scanner  at  an  out  of  focus  loca¬ 
tion.  With  the  system  in  the  scanning  mode,  the  noise  of  the  output 
signal  should  be  less  than  one  quantization  level,  i.e.,  <0.1%. 


3.4.  Response  uniformity 

Response  nonuniformities  can  be  of  two  kinds.  The  response  across 
any  detector  may  be  nonuniform  or,  if  a  detector  array  is  used,  there 
may  be  a  detector-to-detector  and/or  a  gradual  response  variation  along 
the  length  of  the  array. 

The  first  can  be  checked  by  scanning  the  array  across  a  straight 
edge  parallel  to  the  length  of  the  array.  Any  jaggedness  in  the  output 
will  indicate  a  nonuniforra  response,  i.e.,  that  there  is  a  sensitivity 
difference  between  the  leading  and  trailing  edge  of  one  or  more  detectors. 
This  effect  could  introduce  a  significant  mensuration  error  which  would 
not  necessarily  be  apparent  from  a  measurement  of  the  array  geometry. 

A  variation  in  response  along  the  length  of  an  array  can  be  corrected 
in  post-scan  signal  processing;  however,  its  magnitude  has  to  be  known. 

One  method  to  accomplish  this  is  to  scan  the  length  of  an  Eastman  Kodak 
step  tablet,  monitoring  the  output  of  each  detector  at  each  density  level. 
The  two  defects  of  detector-to-detector  variations  or  striping  and  the 
gradual  change  over  several  hundred  detectors  will  be  apparent,  particu¬ 
larly  with  10-bit  quantization.  The  latter  can  be  a  result  of  how  the 
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detectors  were  fabricated  and/or  of  the  uneven  illumination  of  the  target 
to  be  scanned.  Post-scan  signal  processing  should  reduce  these  nonuni¬ 
formities  to  less  than  one  quantization  level,  i.e.,  <0.1%, 


3.5.  Flare  light 

Permissible  flare  light  level  was  not  specified  in  the  AIDS  purchase 
description  and  for  a  well-designed  fixed-aperture  system  may  be  of  no 
consequence.  For  systans  that  provide  various  scanning  aperture  sizes 
and  different  sets  of  apertures,  the  amount  of  flare  or  stray  light  will 
depend  on  the  choice  of  pre-aperture  used  with  the  selected  post-aperture. 
In  badly  designed  and/or  dusty  systems,  the  flare  light  can  be  substan¬ 
tial  and  significantly  reduce  the  MTF.  In  such  cases  it  may  be  possible 
to  enhance  performance  by  introducing  flare  stops  in  the  system. 

The  effects  of  flare  can  be  expressed  quantitatively  in  terms  of  MTF 
reduction  or  arbitrarily  in  terms  of  the  density  measured  when  the  scan¬ 
ning  area  is  blocked  by  a  microdot  2,  5  or  10  times' its  size. 

Figure  23  shows  plots  of  modulation  versus  square  bar  frequency  for 
a  POS  microdensitometer  using  the  Itek  target  described  in  Figure  15. 

The  curves  show  that  increasing  the  size  of  the  pre-aperture  increases 
the  amount  of  flare  light  in  the  optical  system  and  substantially  reduces 
the  MTF. 
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Figure  23:  Scans  to  determine  the  influence  of  stray  light  on  performance 
Conditions:  digital  recording;  lOX  objectives;  white  light; 
200X  magnification;  "A"  apertures;  speed  1  m/s;  1-ym  scan 
aperture.  The  numbers  on  the  graph  indicate  the  pre-aperture 
sizes.  The  dotted  line  is  the  incoherent,  diffraction-limited 
square  wave  MTF  =  500  nm,  NA  =  0.25,  post-aperture  is  1  ym 
square) 
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3.6.  Temporal  stability 

Temporal  stability  Influences  not  only  geometrical  accuracy,  but 
also  radiometric  accuracy.  The  problem  usually  results  from  the  heat 
generated  by  the  light  source  and  electronics  causing  expansion  or  creep, 
during  set-up  and  operation  time,  of  one  part  of  the  optical  system  with 
respect  to  the  other. 

The  PDS  microdensitometer  is,  with  its  pre  and  post-apertures, 
particularly  susceptible  to  thermal  changes  between  parts  of  the  instru¬ 
ment  containing  the  two  apertures.  In  a  12  hour  period  after  turning 
on  the  light  source  the  apertures  moved  apart  by  over  10  pm.  The  change 
in  radiometric  response  is  shown  in  Figure  24  over  a  46  min.  time  period. 
Obviously,  after  the  apertures,  approximately  1  and  2  pm  in  size,  had 
moved  1.5  pm  with  respect  to  each  other,  no  light  was  transmitted  and 
the  density  reading  became  very  high. 


Figure  24:  Density  plot  vs.  warmup  time 
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3.7.  Coherency  effects 

Coherency  effects  are  present  in  most  microdensitometers  and  can 
only  be  minimized  hy  collecting  all  the  diffracted  light  transmitted  by 
the  target.  If  this  is  not  achieved,  density  values  will  be  high  and 
depend  on  the  phase  structure  of  the  target.  Most  processed  film  images 
show  a  relief  image  due  to  changes  in  film  thickness  caused  by  differ¬ 
ential  shrinkage  of  the  image  during  processing. 

Coherency  effects  often  give  rise  to  enhanced  modulation.  Thus  the 
MTF  curve  for  a  perfect  systen  shown  as  the  dotted  line  in  Figure  23  in¬ 
dicates  that  the  2-um  pre-aperature  curve  is  higher  than  the  perfect  curve 
for  spatial  frequencies  below  250  cy/nm.  MTF  curves  generated  this  way 
are  target  dependent  and  only  with  a  phaseless  target  can  the  true  MTF 
be  determined.  Whether  it  is  appropriate  to  measure  the  MTF  with  a  phase¬ 
less  target,  or  with  a  target  produced  on  the  same  film  to  be  routinely 
used  with  the  system,  is  a  matter  for  the  buyer  and  seller  to  argue. 


4.  Test  strategy 

Because  of  the  more  critical  nature  of  the  radiometric  performance 
it  is  suggested  that  the  radiometric  tests  be  performed  prior  to  the 
geometric  tests.  Any  necessary  significant  corrections  to  the  radiometric 
components  should  be  applied  before  entering  the  geometric  test  phase. 

Table  1  gives  a  brief  overview  of  the  radiometric  tests.  Options 
are  indicated  and  priorities  marked. 

The  MTF  determination,  the  check  on  dynamic  range,  response  linearity 
and  uniformity  are  absolutely  necessary.  These  parameters  are  fundamental 
measures  of  the  radiometric  performance  of  the  system.  Flare  light,  tem¬ 
poral  and  coherency  effects  will  influence  those  basic  parameters.  An 
independent  determination  of  flare  light,  temporal  and  coherency  effects 
becomes  important,  only  if  a  thorough  analysis  of  the  system  errors  is  re¬ 
quired. 


45 


The  spatial  resolution  is  defined  by  the  MTF  curve  and  the  specified 
threshold  value  of  0.3  for  the  nrodulus  (AIDS,  1982).  The  resulting 
value  is  based  on  purely  digital  processing  of  the  data.  A  visually 
determined  spatial  resolution  might  significantly  deviate  from  this  value. 

A  threshold  modulation  curve  for  the  display  monitor  would  have  to  be  used 
if  resolution  measurements  are  performed  on  this  monitor. 

Table  2  sunmarizes  the  geometric  tests.  Options  and  priorities  are 
indicated. 

The  global  and  local  gecmietrical  accuracy  tests  are  considered  as 
absolutely  necessary.  If  the  system's  output  is  congruent  to  the  input 
a  3-parameter  transformation  (29)  would  be  sufficient  to  relate  correspond¬ 
ing  input  and  output  features  to  each  other.  Any  additional  significant 
polynomial  term  would  indicate  a  malfunctioning.  For  the  analysis  of 
these  possible  additional  polynomial  terms  the  computation  of  the  poly¬ 
nomial  (1)  for  the  global  test  should  be  done  in  two  different  ways.  One 
computation  should  be  based  on  the  use  of  5  x  5  grid  and  the  true  errors 
should  be  analyzed  at  the  redundant  grid  points.  The  other  should  use 
the  complete  13  x  13  grid  for  polynomial  parameter  determination.  If 
the  polynomial  parameters  of  both  methods  are  identical  (within  the  range 
of  random  effects),  it  is  safe  to  conclude  that  the  system  generates  only 
global  deformations.  Otherwise,  regional  deformations  are  encountered 
and  the  modelling  with  global  polynomials  is  not  optimal. 

A  FORTRAN  subroutine  for  the  computation  of  bivariate  orthogonal 
polynomials,  based  on  a  least  squares  adjustment  as  suggested  in  Chapter 
2.3.1.,  is  available  in  the  Department  of  Geodetic  Science  and  Surveying, 
OSU. 

Another  important  procedure  is  the  check  on  how  fundamental  objects 
like  single  bars,  dots,  and  circles  are  rendered  by  the  system.  This  can 
be  considered  as  a  variation  of  the  standard  resolution  tests.  The  cor¬ 
responding  test  standards  are  selected  such  that  they  can  be  used  with 
variable  scanner  resolution. 


Table  1:  Radiometric  tests 


— 

Test 

purpose 

— 

Test 

standards 

Remarks 

Priority 

1  MTF 

determination 

, 

1 

1 

) 

a)  Itek  15-bar  target 

b)  Razor  blade  edge 
on  Kodak  64 9F 

Measurements  as  function 
of  scan  speed  and  posi¬ 
tion  on  the  format;  re¬ 
cording  in-scan  and 
across-scan  direction; 
for  CCD  linear  arrays: 
test  of  all  detector  out¬ 
puts 

1  Dynamic  range, 

!  linearity, 

I  noise 

I 

a)  Kodak  grey  step 
tabl  et 

b)  Continuous  grey 
scale 

Measurements  as  function  . 

of  scan  speed  and  sampl¬ 
ing  rate 

1 

'  Response 
.  uniformity 

Straight  edge 

Kodak  grey  step 
tablet 

Edge  for  nonuniformity 
across  one  detector; 
step  tablet  for  varia¬ 
tions  along  array  line 

1 

^  1 

Flare  light 

a)  MTF  reduction 

b)  Microdot  block¬ 
ing 

Amount  depends  on  size 
of  pre-aperture 

- 

2 

:  Temporal 
stability 

t  • 

1 

Neutral  low  density 
object 

Compare  density  readings 
over  a  certain  time  peri¬ 
od  after  switching  on 

i  ; 

:  Coherency  | 

effects  1 

1 

a)  Phaseless  tar¬ 
get 

Optimal:  Phaseless  target 

2 

! 

i 

b)  Standard  film 
target 

1 


I 

Table  2:  Geometric  Tests 


Test 

purpose 

Test 

standards 

Remarks 

Wild  standard  test  grid 

a) ,  5x5  grid  for  polynom'^al 

coefficient  detenu! nation; 
check  true  errors  at  other 
grid  points  after  trans¬ 
formation 

b)  13x13  grid  for  polynomial 
coefficient  determination 

c)  3-parameter,  similarity, 
affine  transformation 
with  13x13  grid 

1 

Local 

performance 
and  rendering 
of  fundamental 
objects 

a)  Dot  and  line  com¬ 
parison  chart 

b)  Radial  circle  test 
chart 

Test  standards  designed  for 
variable  scanner  resolution 

j 

1 

:*) 

Skewness  and 
linearity  of 

1  array  line 

i 

a)  Standard  test  grid 

b)  Square  cross-line 
graticule 

Test  lines  inclined  to  scan 
direction;  optimal: 
perpendicular 

2 

i 

*) 

Length  of 
array  line 

1 

i 

» 

Line  test  chart  (no.  32):  Analysis  of  output  of  a 

j  single  array  line 

2 

Stepping  and 
wav i ness 

Standard  test  grid 

Qualitative  check  visually 

1 

2 

1 

1 

Sampi 1 ng  rate  a )  Standard  test  grid 
error  in  scan  .  ,  . 
d1r«:t1or,  ■>> 

! _ 

Check  of  distance  between 
parallel  lines  perpendicular 
to  scan  direction 

2  ! 

1 

i 

1 

Interval 

tracking 

error 

1 

j 

a)  Standard  test  grid  ‘  Test  lines  parallel  to  scan 

b)  Square  cross-line  direction 

graticule  .  Test  chart  no.  32  for  non- 

:  constant  portion 

1 

1 

3 

'  Convergence/ 
divergence 
tracking  error 

Standard  test  grid  Check  of  distance  between 

1  parallel  lines  parallel 
!  to  scan  direction 

3  i 

1 

•  •  • 


applies  to  linear  array  scanners  only 


All  the  other  geometrical  tests  are  largely  designed  for  the  analysis 
of  a  possible  malfunctioning,  and  as  such  for  the  detection  of  specific 
system  errors. 


5.  Conclusions 

For  this  report  the  test  strategy  of  the  AIDS  scanner  had  to  be  set 
up  as  a  device  independent  procedure.  Because  the  final  design  of  the  sys¬ 
tem  was  not  available  all  possible  instrument  errors  could  not  be  antici¬ 
pated.  Consequently  the  tests  had  to  be  constructed  in  a  rather  general 
way.  Therefore  it  was  appropriate  to  develop  options  and  priorities, 
which  are  established  in  Chapter  4. 

It  is  suggested  that  the  radiometric  tests  be  performed  prior  to  the 
geometric  tests.  MTF  determination,  dynamic  range,  noise,  linearity  and 
uniformity  of  the  detector  response  are  considered  as  high  priority  tests. 
The  tests  on  flare  light,  temporal  stability,  and  coherency  effects  are 
second  order  tests.  It  is  crucial  that  the  measurements  for  the  radio- 
metric  tests  are  performed  directly  with  the  digital  data.  The  measure¬ 
ments  for  the  geometrical  tests  should  be  performed  on  the  display  moni¬ 
tor.  High  priority  geometrical  tests  are  those  on  global  performance, 
local  performance  and  rendering  of  fundamental  objects.  The  remaining 
tests  of  Table  2  are  designed  for  the  purpose  of  detecting  specific  in¬ 
strument  errors.  Some  of  these  tests,  including  length,  skewness  and 
linearity  of  array  line,  refer  to  the  application  of  a  linear  array  scan¬ 
ner.  Some  possible  instrument  errors  like  stepping,  waviness,  and  track¬ 
ing  errors  are  easier  to  detect  in  linear  array  systems  than  in  single 
element  type  scanners. 

For  the  implementation  of  these  test  procedures  and  the  processing 
of  data,  some  computer  programs  are  readily  available  from  the  Optical 
Sciences  Center,  University  of  Arizona,  Tucson,  and  from  the  Department 
of  Geodetic  Science  and  Surveying,  The  Ohio  State  University,  Columbus. 
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•3 

The  coordinates  listed  on  the  attached  sheet  are  given  in  millimeters,  rounded  to  10  mm. 

The  coordinate  errors  and  the  mean  square  errors  of  measurement  (computed  from 
differences  of  two  independent  sets  of  measurements)  are  given  in  microns  ( 10  mm), 
rounded  to  the  nearest  tenth  of  a  micron  (10‘^mm). 


GRID  PLATE  KEASURn'EJT 


EDV-NO.:  116765  FLATE-NO.:  7019  DATE:  80/05/13 


FINAL 

COORDINATES  (l-3-l)  AJID 

ERRORS  (NICRONETER) 

POIITT 

X 

DX 

y 

DY 

11 

400.000 

.4 

400.000 

.3 

12 

500.000 

-.5 

400.000 

-.2 

13 

600.000 

-.2 

400.000 

-.4 

21 

400.001 

.6 

499.999 

-.8 

22 

499.999 

-.5 

500.000 

-.5 

23 

600.000 

-.1 

500.000 

.0 

31 

400.000 

.5 

600.000 

.4 

32 

500.000 

-.2 

600.000 

.2 

33 

600.000 

-.0 

600.001 

.9 

I'SE  OF 
MSE  OF 

I'EASURISIT 

CG0RDIII.4TES 

:  NX  « 

:  I'C  = 

.1  MY  = 

.5 

.2 
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Appendix  A. 2:  Square  cross-line  graticules 
(manufactured  by  Heidenhain) 


14 

1  miii^rest>lliM  grid  in  OIAOUR 
OMcripllon:  square  cross-line  grid 
of  100  mm  Side  length,  divided  into 
squares  of  i  mm  side  length. 

Accuracy  of  graduation  ±  15  urn. 

Numerals  are  legible  when  viewed 
through  the  glass. 

Opaque  lines  on  transparent  background. 
GI«m  Mibtlratd:  110X110X3  mm 


15a,  15b 

1  2  mffl-cro«t-lin«  grid  In  OIAOUR 
Oetcription;  squa-e  cross-  ne  gnd 
cl  75  mm  side  .length,  divided  ,nto 
squares  each  of  0.5  mm  side  length. 
Accuracy  of  gradual. on  ±  10  urn 
Numerals  are  legioie  when  viewed 
through  the  glass. 

GI«m  subtlraie:  95  X 85 X  5  mm 


iSa  positive 
opaque  grating  on 
transparent  background 
ISb  negative- 
transparent  grating  on 
opaque  background 


16 

0.1  mm<reas>llne  grid  In  OIAOUR 

Oetcripiion:  square  cross-:. ne  grid 
of  50  mm  side  ength,  divided  into 
squares  each  of  O.l  mm  side  length 
Accuracy  of  graduation  ±  10  urn 
Numerals  are  legible  when  viewed 
through  me  g'ass 
Opaq-e  nes  on  a  f'ansparent 
backgrounq. 

GlaM  Mbilrete:  60X60X5  mm 


jendix  A. 3; 


Dot  and  line  comparison  chart 
(manufactured  by  Heidenhain) 


20a,  20b 


Del  and  line  cemperisen  chart  from 

10  |im  to  1 000  (im.  In  OIAOUR 

as  item  19,  but  with  two  rows  of  19  dots 

and  lines. 

width/dia. 

10,  20, 30, 40,  so,  60, 70, 80, 90, 
100,125,150,200, 250. 

300, 400, 500, 750  and  1000  |im. 

Glass  suhslrate:  60X30X  2  mm 


30a  positive, 
opaque  grating  on 
transparent  background 
30b  negative: 
transparent  grating  on 
opaque  badcground 


lendix  A. 4;  Parallel  line  test  chart 

(manufactured  by  Heidenhain) 


Use  test  chart  In  OIAOUR 

16  groups  with  line  gratings 
from  1.25  imes'mm  up  to  iQO  lines/mm, 
i.  e.  grating  pitoi  from  0.8  mm  up  to 
0.01  mm,  space  to  ime  ratio  ”  i  :  i 
line  length  22  mm 

group  lengths  from  approx.  2.4  mm  up  to 
approx.  5  mm. 

Opaque  lines  on  transparent  background. 

Glass  substrata:  85X35X3  mm. 


'  iW 


5  to  s  7i  « 

X*  >D  S3 


®  I  =  ( 


TraunreuI 
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Appendix  B:  Kodak  high  resolution  aerial  duplicating  film 
CEstar  thick  base)  SO-187 


80>ia7/so-ia2 


KODAK  HIGH  RESOLUTION  AERIAL  DUPLICATING  FILM 
(ESTAB  THICK  BASE)  SO-187 
KODAK  HIGH  RESOLUTION  AERIAL  DUPLICATING  FILM 
{BSTAB  BASE)  SO-192 


Blua-senaltiv*.  axtramaty  fitM^raln  nagativ*  fllma  with  wary  high  raaolving  powar. 
Oaaignad  for  duplicating  fina^aln  aarial  film. 


Baaa: 

SanaMvMy: 

SafailgM: 

Total  nim  Thieknaaa: 

Oal/Baaa  Ratio: 
Watght: 

RMS  Oranularlty  Valua: 
Raaolving  Rowor 

Rodprocity  Ctfact 
AdHMtmofita: 


SO-1S7:  7.0-mil  Estar  Thick  Baaa  with  a  B*11  (fast-drying)  backing. 

SO>1S2:  4.0-mil  Estar  Basa  with  a  B-11  (fast-drying)  backing. 

Blua  sansitiva. 

Usa  a  Kooak  Safalight  Filtar  No.  1A  (light  rad)  in  a  auitabla  lamp  with  a  1 5-watt  bulb  at 
laaat  4  fast  (1.2  matras)  from  tha  film. 

SO-1S7:  Tha  nominal  thieknaaa  of  this  film— including  amulsion  (0.26  mil),  basa 
(7.0  mils),  and  fast-drying  backing  (nil)— is  7.26  mils. 

SO-192:  Tha  nominal  thieknaaa  of  this  film— including  amulsion  (0.26  mil),  basa 
(4.0  mils),  and  fast-drying  backing— is  4.26  mils. 

SO-187:  0.037 
SO-192:  0.06 

80-197:  Tha  weight  of  this  film  (unprocossad)— in  aquilibrium  with  50  percant  relative 
humidity— is  0.054  Ib/lt*  (0.027  g/cm*). 

SO-192:  The  weight  of  this  film  (unprocessed)— in  equilibrium  with  50  percent  relative 
humidity— is  0.032  Ib/ft'  (0.016  g/cm*). 

5  (read  at  a  net  diffuse  density  of  1.0,  46-^m  aperture). 

Test-Object  Contrast  1000:1—800  linas/mm 
Tast-Objact  Contrast  1.6:1 — 250  linas/mm 

For  axposuraa  in  tha  ranga  of  1/1000  to  1/100  sacond,  no  adjustmants  for  raeiprocity 
ara  required. 

For  an  axposure  time  in  the  range  of  1/10  second,  no  adjustment  in  exposure  time  or 
development  time  is  required.  However,  the  lens  aperture  should  be  increased  Vk  stop. 

For  an  expoaure  time  in  the  range  of  1  second,  no  adjustment  in  development  time  is 
required.  However,  either  the  lens  aperture  should  be  increased  1  stop  OR  the  exposure 
time  should  be  increased  to  2  seconds. 


